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Abstract 

The potential V(x) = —x 4 , which is unbounded below on the real 
line, can give rise to a well-posed bound state problem when x is 
taken on a contour in the lower-half complex plane. It is then PT- 
symmetric rather than Hermitian. Nonetheless it has been shown 
numerically to have a real spectrum, and a proof of reality, involving 
the correspondence between ordinary differential equations and inte- 
grable systems, was subsequently constructed for the general class of 
potentials —{ix) N . For such Hamiltonians the natural PT metric is 
not positive definite, but a dynamically-defined positive-definite met- 
ric can be defined, depending on an operator Q. Further, with the 
help of this operator an equivalent Hermitian Hamiltonian h can be 
constructed. This programme has been carried out exactly for a few 
soluble models, and the first few terms of a perturbative expansion 
have been found for the potential m 2 x 2 + igx 3 . However, until now, 
the —x 4 potential has proved intractable. In the present paper we give 
explicit, closed-form expressions for Q and h, which are made possible 
by a particular parametrization of the contour in the complex plane 
on which the problem is defined. This constitutes an explicit proof of 
the reality of the spectrum. The resulting equivalent Hamiltonian has 
a potential with a positive quartic term together with a linear term. 

1 Introduction 

There has been a great deal of interest in non-Hermitian Hamiltonians since 
the numerical observation by Bender and Boettcher £Q that Hamiltonians of 
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the form 

H = p 2 -g(ix) N (1) 

have a real positive spectrum for N > 2. As illustrated in Fig. 1 (from 
Ref. pp, where g — 1), their spectra constitute a smooth extrapolation from 
the simple harmonic oscillator, for which N = 2. The reality of their spectra 
is understood as being due to their unbroken PT symmetry, but there is no 
simple way of telling in advance whether or not this symmetry is broken, as 
indeed it is for N < 2, where the spectra are partly complex. Eventually a 
rather intricate proof of the reality of the spectrum, involving the correspon- 
dence between the differential equations for such potentials and integrable 
models, was constructed by Dorey et al. 
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Figure 1: Energy levels of Eq. (^) for iV = 4 with g = 1, from Ref. jT] 

A potential problem with such Hamiltonians is their physical interpretation, 
since the natural PT norm on the Hilbert space, / dxi/j*(x)ijj(—x), is not 
positive definite, in contrast to the usual norm / dxip*(x)ip(x). However, it 
turns out to be possible to construct an alternative norm, the CPT norm [3], 
which is indeed positive definite. This norm is different from the usual norm, 
in that it is dynamically determined by the Hamiltonian itself, and needs to 
be calculated in each individual case. 
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Such calculations were encompassed by Mostafazadeh [3] in the more general 
framework of pseudo-Hermiticity, whereby 

# f = r]Hr]- 1 (2) 

Here the operator rj is Hermitian and positive definite, and may usefully be 
written as 77 = e~®, in order to connect with the notation of Ref. p2|, where, 
for PT-symmetric Hamiltonians, 77 = PC and Q was defined by C = e®P. 
For calculational purposes it is much easier to deal with Q rather than 77 
directly. Mostafazadeh showed further that 

h = e~5 Q He5 Q (3) 

is an equivalent Hermitian Hamiltonian, obtained from H by a similarity 
(Darboux) transformation. 

In general it is difficult to solve Eqs. (J2J), © exactly; instead one uses per- 
turbation theory in a small parameter e. If H is of the form H = H + eHi, 
where Hq is Hermitian and Hi anti-Hermitian, then Q can be taken as 
Q — Srodd Qr£ r , which then gives h = J2reven h r s r . In this case the first 
few equations for the Q r , arising from the expansion of Eq. read 1 

[Qi,H ] = 2H X (4) 
[Q 3 ,H ] = \[Qi,[Qi,Hi\] 

[Q 5l H ] = 1([Q 3 ,[Qi,Hi]] + [Qi,[Q 3 ,Hi]])-±[Qi,[Qi,[Qi,[Qi,Hi]]\] 

and so on. Using these, the first few equations for the h r , arising from the 
expansion of Eq. 0, can be cast in the form 

^0 = H 

h 2 = -\[Qi,Hi] (5) 
h, = ^[Qi,[Qi,[Qi,Hi]}}-l[Q 3 ,Hi]. 

The smooth continuation from the harmonic oscillator, and the ODE-IM 
correspondence of Ref. rest on the fact that the Schrodinger differential 
equation has several different sectors, defined by wedges in the complex x- 
plane. Along the centre of the wedges the wave-function decays exponentially 
at infinity, while along the edges the wave-function is purely oscillatory. Fig- 
ure 2, taken from Ref. pQ, shows the particular wedge that connects smoothly 
with that for the harmonic oscillator. 

The critical case, where the upper edge of the wedge coincides with the real 
axis is the case N = 4, i.e. the potential — x 4 . For N < 4, it is possible to 
stay on the real axis, where the wave function decays exponentially, albeit 
with an oscillatory modulation, but for iV > 4 we have no option but to 
formulate the problem on a contour in the lower half x plane. 

lr The equations of even order are satisfied identically by Q2n = 0. 
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Figure 2: Wedges in the complex plane in which the Schrodinger equation 
for Eq. (JTJ) is posed. 

This is the fundamental reason why the —gx A potential has proved so in- 
tractable. At first sight it appears Hermitian: it is only because of the 
contour on which it must be defined that it is non-Hermitian. The problem 
is inherently non-perturbative, so any expansion to be attempted can not be 
in the coupling constant g. A previous attempt [Oj used WKB methods, and 
was able to calculate Q to leading non-perturbative order. 

2 Choice of Contour 

Our present approach starts with the idea of Mostafazadeh |7], to map the 
problem back onto the real axis using a real parametrization of a suitable 
contour. A wide variety of contours are possible, as long as they go off to 
infinity at an angle within the wedges. Taking N = 4 and writing the original 
x variable of Eq. (Q) as z to reflect its complex character, so that 



the parametrization used in Ref. [7j was 



corresponding to straight-line contours, with an infinitesimal rounding off 
near the origin. Here 9 was taken as tt/6, the optimal angle for N = 4. The 
resulting Hamiltonian was 



z = x cos 6 — i\x\ sin 9 , 



(7) 
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Because of rounding, there are non-trivial boundary conditions at x = 0, 
namely (i) ip is real and continuous, (ii) ^'(O - ) = e 4z9 ip'(0 + ). Consequently, 
in calculating H* there is an additional term 5H beyond the obvious one. 
A calculation of Q with this Hamiltonian is very difficult because of the 
boundary conditions at x = and the lack of an obvious expansion parame- 
ter. In Ref. jB] we attempted to make an expansion in 9, freeing it from its 
optimal value for N = 4, noting that any positive value for 9 would suffice to 
make the wave function vanish with an exponential component. In addition 
we smoothed out the curve chosen in Ref. [7,, taking the hyperbola 

z = xcos9 — isin9 y/(l + x 2 ), (9) 

in order to remove the boundary conditions at x = 0. Unfortunately this 
calculation did not produce a very useful h, but rather one which still had 
a —x 4 term, so that the asymptotic behaviour of the wave function was 
oscillatory, with only a power suppression. 

In the present paper we adopt a different approach. First we choose a new 
parametrization, whose asymptotes are not in fact in the centre of the wedges, 
but rather are inclined at tt/4 to the real axis, and then we introduce an 
artificial parameter e multiplying Hi, the non-Hermitian part of H(x). 
The contour that turns out to give particularly simple results is of the form 

z=-2iy/(l + ix). (10) 

Notice that with this choice, the PT-symmetry of the original Hamiltonian, 
which is a real function of iz, will be respected by the new Hamiltonian, 
written in terms of x. This new Hamiltonian is in fact 

H = + ix),p 2 } -l p - a (l + i x) 2 , (11) 

where {..,..} denotes the anticommutator, p = d/dx, and for convenience 
we have introduced a = 16g. Separating H into its Hermitian and anti- 
Hermitian parts, and multiplying the latter by the artificial parameter e, 
which at the end will be set equal to one, we write 

H = H + eH 1 , (12) 

where 

H = p 2 - ^p + a(x 2 - 1) 

Hi = - i {x,p 2 } — 2ia x (13) 
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3 Calculation of Q and h 



First we calculate Qi from the first of Eqs. (pEJ), namely [Qx,H ] = 2H\. 
As a general, systematic procedure for such problems we would write the 
Hermitian operator Qi as a sum of anticommutators of the form Qi = 
J2 n0 d([{fn(x),p n }, where f n (x) is a real function of x, and gradually in- 
crease the order n. However, in this case H and Hi are so simple that the 
solution can essentially be found by inspection. Thus a p 3 term in Qi will 
produce the desired structure i {x,p 2 } when commuted with the x 2 term of 
H , while a term in p will produce the x term of Hi. By equating coefficients 
we find that 

3 

Qi = -^ + 2p. (14) 

In order to calculate Q 3 from the second of Eqs. we need the double com- 
mutator [Qi, [Q 1 ,H 1 ]}. First let us calculate the inner commutator [Qi,Hi\, 
which will also be needed for the computation of h: 

4 

[Qi, Hi] = - V — + 4p 2 - 4a. (15) 
a 

The crucial point is that this is a function of p only, and therefore commutes 
with Qi. Thus [Qi, [Qi, Hi}} = 0, which means that Q 3 = 0. Then the third 
of Eqs. (JH) shows that Q5 = and so on. Thus we have an exact solution for 
Q, after setting e — 1, namely 



Q = - ^ + 2p 



(16) 



Having obtained the metric operator Q we are in a position to calculate the 
equivalent Hermitian Hamiltonian h of Eq. (JHJ). Because the expansion for 
Q has truncated, so does that for h, namely h = Hq + h.2- The commutator 
required for the evaluation of h 2 has already been calculated in Eq. (fT5|). so 
it is straightforward to evaluate h, with the remarkably simple result that 



h = Ia~^ P + ax2 ' ^ 

We emphasize that this Hermitian Hamiltonian, defined on the real line, has 
the same energy spectrum as that of the original H of Eq. (jUJ) defined on a 
complex contour. The only unusual feature of h is that it does not have the 
standard form of a quadratic kinetic term plus a potential. However, just 
such a Hamiltonian results if we take the Fourier transform. In terms of the 
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transformed variable y, and after a rescaling y —>■ yy/a, we have 



h = p 2 y + — ay 4 - — Vay 
y 4 2 



4 Discussion 



Equations (fTfij) and (fTSj) constitute our main results. The latter exhibits a 
standard Hermitian Hamiltonian, with a positive quartic potential plus a 
linear term, shown in Fig. 3, whose spectrum is the same as that of the 
original problem, with a —z 4 potential posed on a contour in the complex 
plane. It constitutes the first direct, constructive proof of the reality of the 
spectrum of Eq. (JBJ). In accordance with our introductory remarks, we note 
that h is completely non-perturbative, since, without a harmonic term m 2 x 2 
term in the potential, g can rescaled to 1. 
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Figure 3: The potential of Eq. (JTSJ), with a = 16 (g — 1). 

We have performed a numerical calculation of the energy eigenvalues of 
Eq. (fl~8|) . using both Runge-Kutta integration and the variational truncated 
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matrix method of Ref. j^j. Both methods give eigenvalues that are indistin- 
guishable from those cited by Bender and Boettcher (calculated by Runge- 
Kutta integration along a complex contour) in their original paper pQ. 
A simple extension of the above result can be obtained when an additional 
harmonic term m 2 z 2 is introduced into Eq. (jHJ. The only change in Eq. (|T7jl 
is that h becomes 

h = _1 
4a 2 

with corresponding scaled Fourier transform 

~ h = p y + ^ ay2 ~ 4m ' )2 ~\^ y (20) 

After completion of this work we were made aware of an earlier paper by 
Buslaev and Grecchi which showed the spectral equivalence of the mas- 
sive version of the — x 4 theory (their H e (ig,j), with j = 1), formulated on 
the line z = x — irj, with a Hermitian Hamiltonian that can be identified 
with Eq. (|20|) on setting a = 4g 2 , m = 1/2. Their method made use of 
the perturbation series for the energy eigenvalues of the two Hamiltonians, 
which only exists for m ^ 0. However, they were subsequently able to go 
the massless limit by rescaling and taking g to oo. In this way they obtained 
the spectral equivalence between Eq. © and Eq. (JT7j) (see their Theorem 6, 
with j = 1, a = after a simple rescaling). 

The present paper approaches the problem from a completely different per- 
spective and offers a simple, explicit and transparent derivation of these 
spectral equivalences, together with the operator Q required to define the 
positive-definite metric, and the observables 0] , of the non- Hermitian Hamil- 
tonians. 
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